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Abstract 

A model-independent determination of the primordial power spectrum 
of matter density fluctuations could uniquely probe physics of the very 
early universe, and provide powerful constraints on inflationary models. We 
parametrize the primordial power spectrum A 2 (k) as an arbitrary function, 
and deduce its binned amplitude from the cosmic microwave background 
radiation anisotropy (CMB) measurements of Maxima, Boomerang, and 
DASI. We find that for a flat universe with A 2 s (k) = 1 (scale- invariant) 
for scales k < 0.001 h/Mpc, the primordial power spectrum is marginally 
consistent with a scale-invariant Harrison-Zeldovich spectrum. However, we 
deduce a rise in power compared to a scale-invariant power spectrum for 
0.001 /i/Mpc < k < 0.01/i/Mpc. Our results are consistent with large-scale 
structure data, and seem to suggest that the current observational data allow 
for the possibility of unusual physics in the very early universe. 
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1. Introduction 



The inflationary paradigm is presently the most plausible solution to the problems 
of the standard cosmology Inflation is consistent with all current observational data. 
However, we are still far from establishing a definitive model of inflation. There exists a 
broad range of inflationary models (cf. |Kolb 1997] , [turner 1997] ), many of which appear 
consistent with observational data. In order to quantify what can be known about inflation 
one desires model-independent measurements of the primordial power spectrum of matter 
density fluctuations as these can provide unique and powerful constraints on inflationary 
models. 

Although it has been conventional to take the primordial power spectrum to be 
a featureless power law in analyzing cosmological data, there are both theoretical and 
observational reasons to allow the primordial power spectrum to be a free function. That 
is, some inflationary models predict power spectra that are almost exactly scale-invariant 
( [Linde 19831 ), or are described by a power law with spectral index less than one QFreese 



[Fricman, fc Olinto 1990[ [La fc Steinhardt 1991| ), while others predict power spectra with 
slowly varying spectral indices ([Wang 1994Q , or with broken scale invariance ([Holman et 



|al. 1991aE , Randall, Soljacic, fc Guth 1996 , Adams, Ross, fc Sarkar 1997 , [Lesgourgucsj 
Polarski, fc Starobinsky 1997| , [Lesgourgues 200ti| ). The latter represents unusual physics in 
the very early universe. For example, inflation might occur in multiple stages in effective 
theories with two scalar fields ( [Holman et al. 1991ab| ), or in a succession of short bursts due 
to symmetry breaking during an era of inflation in supergravity models QAdams, Ross, fc 
Sarkar 1997| ). 



There is also tentative observational evidence for a peak in the power spectrum 
of galaxies at k ~ 0.05 Mpc -1 ( [Einasto 1998[ . Baugh fc Gaztanaga 1999 , Rctslaff et al.| 
19981 , IBroadhurst fc Jaffe 2000| , |Gramann fc Suhhonenko 1999) , Gramann fc Hiitsi 2000|) . 



The simplest explanation for such a peak in the galaxy power spectrum is a new feature in 
the primordial power spectrum. 

The cosmic microwave background radiation anisotropies (CMB) are signatures of 
the primordial matter density fluctuations and gravity waves imprinted at the time when 
photons decoupled from matter. The large-scale structure in the distribution of galaxies is 
a direct consequence of the power spectrum of the primordial density fluctuations. Wang, 
Spergel, & Strauss (1999) have explored how one can use the upcoming CMB data from 
the Microwave Anisotropy Probe (MAP; Pennett et al. 19"97| ; |http : //map . gsfc .nasa.gov ) 
and the large-scale structure data from the the Sloan Digital Sky Survey (SDSS; cf., |Gunn 
|fc Weinberg 1996]) to obtain a model-independent measurement of the primordial power 
spectrum, and to extract simultaneously the cosmological parameters. 

In this paper, we implement the concept of a model-independent measurement of the 
primordial power spectrum from Wang et al. (1999) to extract cosmological information 
from the CMB data from Maxima ( [Hanany et al 20"00| ), Boomerang ( |de Bernardis et~al 
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2000), and DASI flHalverson et al. 2001| ). We parametrize the primordial power spectrum 
as a continuous and arbitrary function determined by its amplitude at several wavenumbers 
[which are equally spaced in log(/c)] via linear interpolation. We then measure these 
"binned" amplitudes from the CMB data of Maxima, Boomerang, and DASI. 



2. Measurement of the Primordial Power Spectrum 



We parametrize the primordial power spectrum as 
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K{k) 

A 2 s (k) 
A 2 s (k) 



ki ki_i 



k% ki-i 



ki_i < k < ki, i — 1, n 



k ^ ko — kr. 



where 



kq 



k kri — kr> 



exp 




+ ln(fco) 



0,n 



We impose A 2 s {k) = 1 for k < k = k m i n for two reasons. First, on the largest scales, the 
CMB data is consistent with a scale-invariant primordial power spectrum flSmoot et al.| 
1992| , |Gorski et al. 1996| ), i.e., A 2 s {k) = 1. Second, the bin amplitude of the primordial 



power spectrum on the largest scales is poorly constrained by the CMB data due to the 
effects of cosmic variance. We also assume that A 2 s {k) = a n for k > k n = k max , because k max 
is close to the scale corresponding to the angular resolution of Maxima and Boomerang. 
We choose k min = 0.001 h/Mpc and k max = 0.1 h/Mpc. 

We perform the parameter estimation by computing the CMB angular power spectrum 
Ci(s) for a discrete set of cosmological parameters denoted as s. To save computational 
time and storage space, we have assumed a flat universe and limited our parameter search 
to a grid in six cosmological parameters, {H , Q m , Q b , ai, a 2 , a 3 }, with £7a = 1 — ^m, ^ri = 0. 
The cij (i = 1, 3) parametrize the primordial power spectrum A 2 (k) as in Eq.(fl|) with n = 3. 
Our assumption of a flat universe is consistent with all current observational data, and 
preferred by non-fine-tuned inflationary models. We use the COBE normalization for all 
the theoretical models. 



Following the Maxima, Boomerang, and DASI teams flJafle et al 200 1| , [Lange et al 
2001| , |Pryke et al. 200l|) , we use an offset lognormal likelihood function to define a y 2 



goodness of fit (Pond, Jaffe, fc Knox 2000| ). That is, for a given theoretical model with 
C z (s) =1(1 + 1)C z (s)/(2tt), we define 



X 2 (s) = -2 ln£ 
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In the above equation, D{ is the measured CMB bandpower in the zth bin, Xi is an offset 
which depends upon experimental details, W\ = Bf is the experimental window function 
(where Bi is the beam function). The weight matrix is given by the Fisher matrix 
= — qc q C . ■ The calibration uncertainty is parametrized by u a , the factor of overall 
relative calibration; and which has a dispersion of a a = 0.08 and 0.2 for Maxima-1 and 
Boomerang, respectively. We also consider beam uncertainties for Boomerang. Following 
Lange et al. (2000), we parametrize the beam uncertainty with exp{ — (I + 0.5) 2 [A(# 2 )]}, 
with (A(0 2 )) = 29fwhm AOfwhm — (572. 0)~ 2 , where we have used FWHM = 12.9', and 
A9 FWHM = 1.4'. Hence, we write 



i 
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where we have defined 5 beam = 572.0 2 A(# 2 ), and g(l) = [(I + 0.5)/572.0] ; 



We use the higher precision Maxima and Boomerang data as published in Lee et 



al. 200 1| and [N etterfield et al. 200 1| . The relevant experimental details (beam functions, 
covariance matrix, and log- normal offsets Xj) have not yet been released. It has been shown 
that in the absence of the knowledge of Xi, Gaussian statistics describe the data better 



than a log- normal distribution with arbitrary Xi (|Bond, Jaft'e, fc Knox 2000]) . Therefore, we 
replace Eq.(H) with 
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where C l data are the experimental band-powers with measurement errors <Tj, and 
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where Ofwhm — 10' and 12.9' for Maxima and Boomerang respectively. Eq. (|T0|) assumes 
symmetric error bars on C data . For the asymmetric errors given by Maxima, we take 
(Tj = <Tj _ for C l data > C l BP (s), and a, = a ij+ otherwise. 

The DASI team has released the experimental details (window functions 
fu, log- normal offsets Xi, and covariance matrix Vy) together with their data 
( |http : / /astro . uchicago . edu/ dasi/| ) . This enables us to use log-normal statistics [Eqs. (|J)- 
(|)], with the weight matrix given by the inverse of the matrix = + a a DiDj, 
where a a = 0.08 is the calibration uncertainty of DASI. 

For a given set of cosmological parameters, {H , Q m , Qb, a±, a 3 }, we marginalize 
the model over the calibration uncertainties of Maxima, Boomerang, and DASI, and over 
the beam uncertainty of Boomerang. The beam uncertainty is only marginalized for 
Boomerang, since only this experiment appears to have significant beam uncertainty. 



3. Results 

Since only five values of H = 100/ikm/sMpc- 1 (h = 0.5, 0.6, 0.7, 0.8, 0.9) are 
computed for our grid, we present our results for the different values of h separately, instead 
of marginalizing over h. 

Figs.l(a)-(c) show the data from (a) Maxima, (b) Boomerang, and (c) DASI, together 
with best-fit models. In each figure, the solid curve is the best-fit model to the combined 
Maxima, Boomerang, and DASI data, while the dotted line is the best-fit model to the 
(a) Maxima, (b) Boomerang, and (c) DASI data separately. The dotted error bars are the 
errors on each data point including calibration uncertainty. We have combined the data 
from Maxima, Boomerang, and DASI assuming that the three experiments are independent 
of each other. 

Table 1 shows the best-fit models to the combined Maxima, Boomerang, and DASI 
data. 

Table 1: Best fit models to the combined Maxima, Boomerang, and DASI data 
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0.9 
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0.05 
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0.9 
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0.990 


0.05 


29.0687 



- 6- 



Table 2-4 show the best-fit models to the Maxima, Boomerang, and DASI data 
separately. 



Table 2: Best fit models to the Maxima data 



h Vt m Qb Qj. Q2 Q3 U mx Xmin 

0.5 0.800 0.100 3.0 1.5 1.5 0.997 4.3313 

0.6 0.600 0.100 2.5 1.0 1.5 0.997 3.2590 

0.7 0.400 0.080 0.9 0.6 1.0 1.010 3.9188 

0.8 0.300 0.070 1.0 0.6 1.2 0.994 4.5866 

0.9 0.200 0.060 0.7 0.5 1.3 1.000 4.7581 



Table 3: Best fit models to the Boomerang data 



h Vt m Vtb Oj. Q-2 Q3 ^£00 dbeam Xmin 

0.5 0.600 0.050 5.0 3.0 1.5 0.968 -0.05 11.4684 

0.6 0.375 0.050 2.5 1.5 1.0 0.968 0.00 11.7898 

0.7 0.250 0.040 2.0 1.2 0.9 0.984 0.00 12.1868 

0.8 0.150 0.030 1.5 0.9 0.8 1.016 0.00 12.4757 

0.9 0.120 0.030 0.9 0.8 0.9 1.008 0.05 13.5643 



Table 4: Best fit models to the DASI data 



h Vt m i\ Qj a 2 Q3 UDASI Xmin 

0.5 0.800 0.070 2.5 1.7 1.3 0.990 4.3187 

0.6 0.450 0.050 4.0 1.7 1.5 1.010 3.3240 

0.7 0.275 0.040 3.0 1.2 1.3 1.006 2.5608 

0.8 0.200 0.040 2.0 0.9 1.3 0.994 2.6985 

0.9 0.120 0.030 0.9 0.6 1.0 1.029 2.3055 



We constrain parameters individually by marginalizing over all other parameters. 
Figs. 2 (a)- (e) show the likelihood functions for the combined Maxima, Boomerang, and 
DASI data for the set of parameters {fi m , Q b , a±, 02, ^3}, for h = 0.5 (dot), 0.6 (solid), 0.7 
(dashed), 0.8 (long-dashed), and 0.9 (dot-dashed). The likelihood functions have been 
derived using \ 2 values which resulted from multi-dimensional interpolations of the \ 2 
values computed for the grid of models (cf. |Tegmark fc Zaldarriaga 2000a|) . 
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Fig.3 shows the primordial power spectrum A 2 (k) measured from the combined 
Maxima, Boomerang, and DASI data for h = 0.6 (solid), and 0.7 (dotted). The ±lcr errors 
are deduced from Fig.2(c)-(e). 

Note that A 2 (k) = 1 corresponds to the scale-invariant Harrison-Zeldovich 
n s — 1 spectrum, with the primordial scalar power spectrum conventionally defined as 
k A 2 (k) oc k na . Clearly, the primordial power spectrum is marginally consistent with a 
scale-invariant Harrison-Zeldovich spectrum. However, there is a rise in the power at 
0.001 /iMpc -1 < k < 0.01 /iMpc -1 compared to a scale-invariant power spectrum. Although 
this rise is not statistically convincing, it could be an indication of interesting physics 
worthy of future investigation. 



4. Discussion 



Due to parameter degeneracies between h and Q m and Qb, the CMB data constrain h 
only weakly. Fortunately, there are a number of independent methods for constraining h 
( [Freedman et al 200 1| , [Branch 1998| , |Kundic et al 1997] ). The current combined CMB data 
from Maxima, Boomerang, and DASI marginally prefer h = 0.6 (see Table 1). 

Although the CMB data is sensitive to Q m h 2 and Q^h 2 , the likelihood curves for 
different values of h do not overlap [see Figs. 2 (a) (b)], since these curves correspond to 
different values of I^a (we assume that fl^ = 1 — f2 m ). The likelihood curves for Qt,h 2 are 
less sensitive to the value of h [see Fig. 2(b)]. This is because varying Qb has a much smaller 
effect on Q\ than varying Q m . 

Our measurement of the primordial power spectrum [see Fig. 3] is consistent with the 
large-scale structure data which seem to indicate a peak in the matter power spectrum at 
k ~ 0.05 Mpc^ 1 (lEinasto 1998| , Baugh fc Caztanaga 19991 , |Retslafi et al. 1998| , Broadhurst 
& Jaffe 



Cramann & Suhhonenko 19991, ICramann & Hiitsi 



We note that the 



real space power spectrum of the PSCz Survey, from 0.01 to 300 h/Mpc ( Hamilton 
l'egmark 2000|), does not show this peak. The data from 2df (jDalton et al. 2000|) and 



SDSS (|Cunn fc Weinberg 1996| ) should help clarify any features in the large-scale structure 
power spectrum. However, if there is a feature at k < 0.002 /iMpc -1 (corresponding to the 
characteristic length scale of the SDSS, see |Wang, Spcrgcl,fc Strauss 1999|) , satelite CMB 
data from MAP or Planck will be required to constrain such features in the primordial 
power spectrum ( |Wang, Spergel,fc Strauss 1999| ). 



A number of authors have used Maxima and Boomerang data to derive cosmological 



constraints (Abazajian, Fuller, & Patel 2000, Amendola 2001 


Avelino et al. 2000, Balbi et 


al 2000, Bento, Bertolami, & Silva 2001 


, Bouchet et al. 2002 


. Brax, Martin, & Riazuelo 


2000, Bridle et al. 2001, Contaldi 200C, Durrer & Novosyadlyj 2001, Enqvist, Kurki-Suonio, 


k Valiviita 200q, Fsposito et al. 2001 


Griffiths, Silk, & Zaroubi 2001, Hannestad 



2000 


, biannestad & Scherrer 2001 


, Hu et al. 2001, Jaffe et al 2001, Kanazawa et al. 


2000 


, Kinney, Melchiorri, Riotto 2001, Landau, Harari, & Zaldarriaga 2001, Lange 


et aJ 


. 2001, Lesgourgues & Peloso 2000, McGaugh 200q, Melchiorri & Griffiths 


2001 


, Padmanabhan & Sethi 2001 


, Tcgmark & Zaldarriaga 2000b, Iegmark, Zaldarriaga, 


& Hamilton 2001, White, Scott, & Pierpaoli 2001|). Our work is unique in allowing the 



primordial power spectrum to be an arbitrary function, thus allowing the possibility for 
detecting new features in the primordial power spectrum. 

Our results seem to indicate that the current observational data do not rule out 
unusual physics (such as multiple-stage inflation) in the very early universe. The upcoming 
data from the CMB satellite missions MAP flBennett et al. 1997| ) and Planck ( De Zotti et 



al. 1999|) , and the large-scale structure data from 2df (Palton et al. 2000|) and SDSS (|Gunn 



fc Weinberg 1996]) should allow for a more definitive measurement of the primordial power 
spectrum (|Wang, Spergel,fc Strauss 1999] ) . These data will more precisely constrain the 
possibility for such complex physics in the very early universe. 
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Fig. 1. — The data from (a) Maxima, (b) Boomerang, and (c) DASI, together with best- 
fit models. In each figure, the solid curve is the best-fit model to the combined Maxima, 
Boomerang, and DASI data, while the dotted line is the best-fit model to the (a) Maxima, 
(b) Boomerang, and (c) DASI data separately. The dotted error bars are the errors on each 
data point including calibration uncertainty. 

Fig. 2. — The likelihood functions for the combined Maxima, Boomerang, and DASI data 
for the set of parameters {fl m , ^b, Q-i, 02, 03}, for h = 0.5 (dot), 0.6 (solid), 0.7 (dashed), and 
0.8 (long-dashed). 

Fig. 3. — The primordial power spectrum A 2 s (k) measured from the combined Maxima, 
Boomerang, and DASI data for h = 0.6 (solid), and 0.7 (dotted). The ±la errors are 
estimated from Fig.2(c)-(e). 
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Fig. 1. — The data from (a) Maxima, (b) Boomerang, and (c) DASI, together with best- 
fit models. In each figure, the solid curve is the best-fit model to the combined Maxima, 
Boomerang, and DASI data, while the dotted line is the best-fit model to the (a) Maxima, 
(b) Boomerang, and (c) DASI data separately. The dotted error bars are the errors on each 
data point including calibration uncertainty. 
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Fig. 1. — (b) Boomerang data. 
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Fig. 1.— (c) DASI data. 




Fig. 2. — The likelihood functions for the combined Maxima, Boomerang, and DASI data 
for the set of parameters {fl m , fib, ai, a 2 , a 3 }, for h = 0.5 (dot), 0.6 (solid), 0.7 (dashed), and 
0.8 (long-dashed), (a) Vt m h 2 . 
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Fig. 2.- (b) tt b h 2 . 




Fig. 2.— (c) ai. 
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Fig. 3. — The primordial power spectrum A 2 s (k) measured from the combined Maxima, 
Boomerang, and DASI data for h = 0.6 (solid), and 0.7 (dotted). The ±la errors are 
estimated from Fig.2(c)-(e). 



